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Abstract 

We develop a systematic framework for exclusive rare B decays of the type B — )• 
K^*h^l^ at large dilepton invariant mass q^. It is based on an operator product 
expansion (OPE) for the required matrix elements of the nonleptonic weak Hamil- 
tonian in this kinematic regime. Our treatment differs from previous work by a 
simplified operator basis, the explicit calculation of matrix elements of subleading 
operators, and by a quantitative estimate of duality violation. The latter point is 
discussed in detail, including the connection with the existence of an OPE and an 
illustration within a simple toy model. 
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1 Introduction 



The rare decays B K'^*H+l- are among the most important probes of flavour physics. 
They are potentially sensitive to dynamics beyond the Standard Model (SM) and have 
been intensely studied in the literature [1]. Measurements have been performed at the 
5-meson factories [2H6] and at the Fermilab Tevatron [7] . Excellent future prospects for 
detailed measurements are provided by the LHC experiments ATLAS, CMS, and LHCb 
at CERN [1], and, in the longer run, by Super Flavour Factories based on e^e~ colliders 
[HHn]. 

The calculability of -B — > K^*h'^l~ decay rates and distributions benefits from the fact 
that these processes are, to first approximation, semileptonic modes. Correspondingly, 
the hadronic physics is described by i? — > K^*'^ form factors, which multiply a perturba- 
tively calculable amplitude. This simple picture is not exact because also the nonleptonic 
weak Hamiltonian at scale ~ has B — >■ K^*U'^l~ matrix elements. The prominent 
example is given by hadronic interactions of the form (s6)(cc), where the charm quarks 
annihilate into through a virtual photon. Such charm-loop contributions are more 
complicated theoretically than the form-factor terms. Even though the charm loops are 
subdominant numerically in the kinematical regions of interest, they cannot be com- 
pletely neglected. In particular, the related uncertainty needs to be properly estimated 
in order to obtain accurate predictions. 

We need to distinguish three regions in the dilepton invariant mass g^, for which 
the properties of charm loops are markedly different. For 7GeV^ ^ q'^ ^ 15 GeV^ the 
presence of very narrow cc resonances leads to huge violations of quark hadron duality 
[12] and the hadronic backgrounds from B — )• K^*^ip, followed hj ip ^ dominate 
the short distance rate by two orders of magnitude. This region in can be removed 
by experimental cuts. 

For < 7GeV^ the kaon is very energetic and the charm loops can be computed 
systematically in the heavy-quark limit using QCD factorization for B decays into light- 
like mesons [TSUI]. This approach was first employed for B — )■ K^*H'^l~ in [TJ]. The 
results have many applications. A summary with detailed references can be found in pQ 
(see also |16j for a recent analysis). 

The high-g^ region, g^ > 15 GeV^ has received comparatively little attention. In this 
case the kaon energy is around a GeV or below, and (soft-collinear) QCD factorization is 
less justified, becoming invalid close to the endpoint of the spectrum at g^ = (m^ — m^^)^. 
On the other hand, the large value of g^ defines a hard scale for the hadronic contri- 
bution to i? — K^*H^l~. Consequently an operator product expansion (OPE) can be 
constructed, which generates an expansion of the amplitude in powers of EK/^/(f (or 
^qcd/ \f<f)- Charm loops, and other hadronic contributions, are thus approximated as 
effective interactions that are local on the soft scales set by and A.qcd- This sim- 
plifies the computation substantially. In fact, to leading order in the OPE the hadronic 
contribution reduces to a standard form-factor term. This picture has been first dis- 
cussed at lowest order in the OPE in [I7j, where it was applied to the endpoint region 
of 5 ^ K^*h+r and B Knl+r. In [H] the OPE was considered in some detail, 
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including a discussion of power corrections. 

In the present paper we formulate the OPE for the high-g^ region of 5 ^ K^*H+l- 
from the outset. Although our approach is similar in spirit to the analysis of [18\, the 
concrete implementation is different. We will also go beyond the estimates presented 
in [18] in several ways. An important difference is that combines the OPE with 
heavy-quark effective theory (HQET), whereas we prefer to work with 6-quark fields 
in full QCD. The latter formulation has the advantage of a simplified operator basis, 
which makes the structure of power corrections and their evaluation considerably more 
transparent. We also retain the kinematical dependence on in the coefficient functions, 
rather than expanding it around = m^. We further discuss the issue of quark-hadron 
duality, which appears relevant because of the existence of cc resonance structure in 
the region of interest. Violations of duality are effects beyond any finite order in 
the OPE. Using a resonance model based on a proposal by Shifman, we quantify for 
the first time the size of duality violations in the high-g^ region of i? — )■ K^*H'^l~ . 
Further aspects and new results of our analysis will be summarized in sections 8 and 
9. The main conclusion is that B — )■ K^*H~^l~ is under very good theoretical control 
also for g^ > 15GeV^. Precise predictions can be obtained in terms of the standard 
form factors, with essentially negligible effects from the additional hadronic parameters 
related to power corrections and duality violation. 

The paper is organized as follows. Section 2 collects basic expressions for later refer- 
ence. In section 3 our formulation of the OPE for B K^*Hn- at high g^ is described 
and the power expansion is constructed explicitly, complete to second order in 1/ ^/q^ 
and with a discussion of weak annihilation as an example of a (small) third-order correc- 
tion. In section 4 we present an estimate of the matrix elements of higher- dimensional 
operators and quantify their impact on the decay amplitudes for both B K and 
B — )■ K* transitions. Section 5 discusses the connection between the OPE for large g^ 
and QCD factorization for energetic kaons, which are shown to give consistent results 
at intermediate g^ ~ 15GeV^. In section 6 we address the subject of duality violation 
in the context of a toy model analysis. The estimate is then adapted to the case of 
B —7- Kl^l" in section 7. In this section we also address conceptual aspects relevant for 
the existence of the OPE and the notion of quark-hadron duality. A comparison of our 
approach with the literature is given in section 8 before we conclude in section 9. Details 
on the basis of operators in the OPE are described in appendix A and some numerical 
input is collected in appendix B. 



2 Basic formulas 



2.1 Weak Hamiltonian 

The effective Hamiltonian for h — )■ sl^l~ transitions reads fT9]|20ll21 

G 



n 



eff 



p=u,c 



CiQi + C2Q2 + CiQi 



i=3,...,10 



(1) 
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where 

K = v;,v,, (2) 

The operators are given by 

i{sp)v-A , Ql = {Pibj)v^A{SjPi)v-A , 

1 {m)v-A , Qi = {Sibj)v-A {qjqi)v-A , 

1 {m)v+A , Qe = {sihj)v-A 



Ql = 


(pb) 


Q3 = 


(sb) 


Q,= 


(sb) 


Q7 = 


e 






87r2 


Q9 = 


a 


2^^ 



i)V+A , 



(3) 



OTT 



a _ - 

v-a{U)v, Qio = ^{sb)v-A{ll)A 

/vr 

Note that the numbering of reversed with respect to the convention of [19]. Our 

coefficients Cg^io correspond to Cg^io in [12] and we include the factor of a/{2n) in the 
definition of Qg^io- The sign conventions for the electromagnetic and strong couplings 
correspond to the covariant derivative = d ^+ieQ f A^+igT"" A'^^. With these definitions 
the coefficients Cj^^ are negative in the Standard Model. 

2.2 Dilepton-mass spectra and short-distance coefficients 

We define the kinematic quantities s = q'^/m^B (where q^ is the dilepton invariant mass 
squared), rx = m\/m\, and 

\K{,s) = l+r\ + s^ ~2rK-2s-2rKS (4) 

The differential branching fractions for B — )■ Kl^l^ can then be written as [22] 

^/^"^'^ = rs^f^WM^ . Xf{s)flis) {\a,iKll)\^ + MKlir) (5) 

The coefficient ag{Kll) contains the Wilson coefficient Cg{fi) combined with the short- 
distance parts of the B — )■ Kl^l~ matrix elements of operators Qi,...,Q8- The co- 
efficient ag{Kll) multiplies the matrix element of the local operator Qg in the decay 
amplitude. The coefficient aio{Kll) = Cio of the operator Qiq is determined by very 
short distances ~ 1/Mw and is precisely known. 

The corresponding formulas for B — K*l^l~ can for instance be found in |23j . 



3 OPE for B Ml+l- amplitudes at high 
3.1 General structure 

The amplitudes for the exclusive decays B — Ml^l" , where M = K, K*, or a similar 
meson, are given by the matrix element of the effective Hamiltonian in ([T]) between the 
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initial B meson and the Ml^l~ final state. The dominant contribution comes from the 
semileptonic operators Qg^io- Their matrix elements are simple in the sense that all 
hadronic physics is described by a set of -B — > M transition form factors. This is also 
true for the electromagnetic operator Qj. The matrix elements of the hadronic operators 
Qi, • • • , Qe, Qs are more complicated. They are induced by photon exchange and can be 
expressed through the matrix element of a correlator between the hadronic part of the 
effective Hamiltonian 

6,8 

HP = C^Q\ + C2Ql + Y,C,Q^ (6) 

i=3 

and the electromagnetic current of the quarks 

f^Q^q^y (7) 

where Qq is the electric charge quantum number of quark flavour q and a summation 
over q is understood. The decay amplitude may thus be written as 



AiB MI+1-) = -^—Xt 



(8) 



where u and v are the lepton spinors and 



Q 2 /* 

A>'^=C,{M\-sY{l-^,)h\B)-^i J d^xe'"-^ {M\T f{x)H%0)\B) 

+ Cr^qx {M\sa"'{l + ^,)b\B) 
q^ 

"^cu = ^'l {M\Tf{x){H\0) - H%0))\B) 

A'l, = Cio{M\sYa-l5)b\B) (9) 

For 6 — s transitions the contribution from y4(^^ is suppressed by the prefactor Xu/Xt and 
can be neglected. 

Exploiting the presence of the large scale ~ ml, an operator product expansion 
(OPE) can be performed for the non-local term 

iCiriq) = — / d^xe*'?-^Tj'^(a;)if'=(0) (10) 
J 

which describes the contribution of 4-quark operators to the b — )■ sl~^l~ amplitude. 
Such an OPE corresponds to integrating out the hard quark loop, leading to a series 
of local effective interactions for the high-g^ region. To leading order in the large-g^ 
expansion this has been presented in [T7j. A discussion of the OPE including higher- 
order contributions has been given in |18j . 
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Figure 1: OPE for the correlation function /C^: Leading-power contributions (operators 
of dimension 3). The sohd square and the virtual-photon attachment indicate the in- 
sertion of the weak Hamiltonian and of the electromagnetic current, respectively. The 
diagram on the left shows the lowest-order term in QCD. On the right is a sample 
diagram for the next-to-leading QCD corrections of order a^. 

Before going into more detail we discuss the basic structure of the OPE for /C^. The 
expansion may be written as 



The operators Od,n are composed of quark and gluon fields and have the flavour quantum 
numbers of [sb). They are ordered according to their dimension d and carry an index n 
labeling different operators with the same dimension. The Cd,n{(l) are the corresponding 
Wilson coefficients, which can be computed in perturbation theory. The large scales 
justifying the expansion are and q^. They are counted as quantities of the same 
order. The coefficients then scale as Cd,n ~ in the heavy-quark limit. Since the 

matrix elements {M\Od,n\B) scale as ^/rr^, the matrix element of each term in (JTTll 
behaves as m^'^~'^ . Current conservation implies that all operators are transverse in g. 



It is convenient to work with the 6-quark field in full QCD. This field could be 
further expanded within heavy-quark effective theory (HQET), in order to make the m^- 
dependence fully explicit. In such an approach many additional operators would arise 
whose hadronic matrix elements are not readily known. In contrast, the advantage of 
using the 6-field in full QCD is that fewer operators appear and that the matrix elements 
of the leading ones are given by common form factors. In this method the OPE becomes 
particularly transparent and we will adopt it here. 

At leading order in the OPE (c? = 3), illustrated in Fig. [H and in the chiral limit 
[rus = 0) there are two operators 




(11) 



d,n 



%oi^ ^ 



(12) 




(13) 




(14) 



5 



Figure 2: OPE for the correlation function /C^: Second-order power corrections (opera- 
tors of dimension 5). The crossed circles denote the places where the virtual photon can 
be attached. 



Using the equations of motion for the external quarks it can be shown that all possible 
bilinears s^Ffe and arising from the correlator /C^ can be expressed in terms of 

(fT3|) and (fT4|) . Consequently, no independent dimension-4 operators of the form 

can appear in the OPE. The complete proof is given in appendix A. As an example, the 

operator si^^{l + 75)6 satisfies the equations-of-motion identity (for = 0) 

M^'^il + 75)& = -^s^il - l,)b + Id^isa^^il + l,)b) + ^9^(s(l + 75)^) (15) 
For any B ^ Xg matrix element with momentum transfer q this is equivalent to 

s^^^(l + 75)6 = -^S7^(l - 75)& - '^qusa'-'il + l,)b + \qn{l + 75)^ (16) 

Because of current conservation only the transverse part {g^\ — q^qx/q^)0^ of such an 
operator can appear in the OPE. From (fT6|) we see that this part can be reduced to 
a linear combination of (fT3l) and ( fT4l) . 

If we keep rris 7^ 0, two additional operators have to be considered 

Oi, = (^g'^'' - ^7,(1 + 75)& (17) 

Ol, = '^q,-sa"^{l~j,)b (18) 

Since ms/mi, is small, and numerically similar to A/rUb, we may formally count these 
as operators of dimension 4. Because they are absent at order their impact will 
be suppressed to the level of asms/rrib ~ 0.5%, which is negligible. Note that these 
operators do in any case not introduce new hadronic form factors. 

At d = 5 (Fig. [2]) we encounter operators with a factor of the gluon field strength 
G^u, which have the form 

= sigGT^rb (19) 

where the F^ denote Dirac and Lorentz structures. We will treat the OPE explicitly to 
the level of d = 5, that is including power corrections up to second order in A/rrib. 
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Figure 3: OPE for the correlation function /C^: Weak annihilation as an example of 
third-order power corrections (operators of dimension 6). Crossed circles denote the 
various virtual photon attachments. 



Although we will not give a full treatment of dimension-6 corrections, we consider as 
an example the effect of weak annihilation (Fig. [3]). This contribution is characterized 
by the annihilation of the two valence quarks in the B meson in the B ^ M transi- 
tion through the weak Hamiltonian. It is described by 4-quark operators, which read 
schematically 

(^6^ann,n = (^ri6sr2r)(^ (20) 

with Lorentz and Dirac structures indicated by (Fi 2)^, and the light quark field r = u, 
d in the case of non-strange B mesons. Weak annihilation provides a mechanism to 
break isospin symmetry, directly at the level of the transition operator. In /C^ weak 
annihilation, in addition to being a third order power correction, comes only from QCD 
penguin operators, which have small coefficients. The contribution to isospin breaking 
from this source will therefore be strongly suppressed. 



3.2 OPE to leading order in CKg 

In this section we give explicitly the first few terms in the OPE to leading order in 
renormalizat ion-group improved perturbation theory, that is neglecting relative correc- 
tions of 0{as)- This order for /Cj^^ is relevant in the next-to-leading logarithmic approx- 
imation to the B — )■ Ml~^l~ amplitude. We may then write 



(21) 



The lower indices of the terms on the r.h.s. denote the dimension d of the corresponding 
local operators, which come with a coefficient of order l/m^~^. The first term reads 



S7^(l -75)&' 



h{z, s){Ci + 3C2 + 3C3 + C4 + 3C5 + Ce) - -h{l, s)(4C3 + 4C4 + 3C5 + Cg) 



-h{0, s){Cs + 3C4) + ^(3^3 + C4 + 3^5 + Ce) 



(22) 
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The coefficient in (122!) requires a UV renormalization, which has to be consistent with 
the definition of Cg. The expression given here corresponds to the NDR scheme used in 
|19j . The function h{z,s) is {z = mc/mh, s = q'^/ml, x = Az'^/s = 4m^/g^) 



h(z, s) = — In — - - In z H h -a; + -(2 + x)^/l - x In \ + zvr (23) 

^ ^ 9 /i 9 27 9 9^ ' V l + VT^ ' 



Next we have 



Here 



X 



1 1- v^l^ , • 
m , + iTX 



y/l^ V 1 + 



X 



(24) 



(25) 



with X = 4m^/g^. The charm-loop contribution in (12^ . proportional to Ci, can be 
inferred from [23]. Note that here we use the convention e^^^s _ _-|^^ j^-^ writing we 
have neglected terms with the small QCD penguin coefficients C^, . . . ,Cq. 
Finally, weak-annihilation diagrams give the dimension-6 term 



r=u,d 



2Qr{fa''il - lb)hj Skl^il - 75)n ~{^^^ A}) 



-ie^'^^^'fi-ipil - 75)&, Sfc7,(l - 75)n 



+ 



3 



E 

r=u,d 



Qr{fi{l - l5)bj SkCT^^l + 75)n + fia^^l - 75)6j Sk{l + 75)n) 



-^(f,(l - 75)6j SfcCT^^l + 75)n - r^cT^^l - 75)&i Sfc(l + 75)^) 



(26) 



The terms in fl26|) only arise from QCD penguin operators, which have small coefficients. 

We remark that all operators in ( !22l) . ( l24l) and ( l26|) vanish identically when contracted 
with g^, as required by gauge invariance. 



3.3 0{cx.s) corrections to the charm loop 

The non-factorizable 0{as) corrections to the charm loop arise from diagrams like the 
one shown on the r.h.s. of Fig. [H The g^-dependence has been recently calculated in 
analytic form as a Taylor expansion in the small parameter z = ml/ ml [25]. Analytic 
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results for rric = had been presented in [26]. In the kinematical range relevant to 
our considerations, it has been shown that the convergence of the series is very good. 
We therefore use the mathematica input files provided by the authors of [25] in the 
online preprint publication for a numerical estimate. We find that the non-factorizable 
0{as) corrections to the charm loop lead to a 10-15% reduction of the real part of ag 
and contribute a negative imaginary part of again 10-15% relative to the short-distance 
contribution from Cg (the precise value is scheme- dependent). This is in agreement with 
the effect found for the inclusive B — )■ Xgl^l" rate in the high-g^ region, as discussed 
in [25], and is similar to the effect observed for the low-g^ region in the exclusive decay 
modes, see Table 5 in [27] . 

It is to be stressed that these corrections almost compensate the factorizable charm- 
loop contribution (diagram on the l.h.s. in Fig. [1]). The reason why the 0{as) corrections 
are not suppressed stems from the different colour structure of the diagrams. Whereas 
the factorizable charm loop comes with a colour-suppressed combination of Wilson co- 
efficients, the additional gluon exchange allows the cc-pair to be in a colour-octet state 
with no such suppression. At even higher orders in perturbation theory, 0(a") with 
> 2, on the other hand, the numerical effect on ag should really be small, as no new 
additionally enhanced colour structures will arise. 

4 Matrix elements and power corrections 

The computation of the amplitude from the OPE requires the evaluation of the matrix 
elements of the local operators. We estimate in particular the matrix elements of the 
higher-dimensional contributions. This will allow us to quantify power corrections to the 
B — ?► K^*^l^l~ amplitude at high g^. The cases of B ^ K and B — t- K* transitions will 
be considered in turn. 



4.1 B ^ K 



The matrix element of the leading dimension-3 operator is given in terms of the familiar 
form factors f±, defined hy {p = k + q) 

{K{k)\-sr{l - l,)h\B{p)) = 2U{q') P + [U{q') + Uq')] q^ (27) 

At the level of the dimension-5 correction in (12^ one encounters operators of the 
form sGap^h. Their matrix elements introduce, in general, new nonperturbative form 
factors. Using Lorentz invariance and the antisymmetry of G^/s and a'''^ one can show 
that 

q^ {K{k)\-sG^p{g^'^a^^ - g'^^a^'^){l + i^)h\B{p)) = (28) 

In order to estimate the remaining term we assume A <^ Ek ttlb for the kaon energy 
Ek- In this limit the matrix element can be computed in QCD factorization. To leading 
order we then find 

k ■ q 



{K{k)\JCUBip)) 



-CiQcj{x)- 



N 



(29) 
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where Cp = {N"^ — 1)/(2A^), iV is the number of colours, and l/A^ is the first inverse 
moment of the B-meson hght-cone distribution amphtude. This matrix element scales 
as (A/mfo)^ relative to ( l27l) in the heavy quark limit, and as 0{\) for large A^. 
In a similar way we can estimate the weak annihilation term 



qZ 



(30) 



This contribution is power-suppressed as {K/mi,Y relative to f l27|) . The suppression by 
the small Wilson coefficients 6*3^4 is partly compensated by a large numerical factor of tt^. 
To relative order {A/mb)^ there is no contribution from the term with C5 and Cq. Note 
that the result in (l30ll also corresponds to the matrix element obtained when naively 
factorizing the four-quark operators. 

Normalized to the amplitude coefficient ag = Cg + . . ., the power corrections from 
)Ch5 and JCnea read 

Aag,,,{K) = ^]^C,QJ{x)j-j-^^ (31) 

Aag,,UK) = -{c, + ^)'-^j^^ (32) 

where r = u, d refers to the spectator quark in the B meson. 

Numerically, we find Aag^HbiK) = 0.019 — 0.012? at = 15 GeV^ for central values 
of the parameters. This number comes with a substantial uncertainty, in particular 
from Xb and /+. Nevertheless, the correction to ag ~ 4 is very small, most likely 
below 1% in magnitude. The correction in ( 13T]) diminishes further for larger g^, reaching 
Aag^H^iK) = 0.006 — 0.002z at the endpoint. Towards the endpoint the kaon becomes 
soft and the result in (!3T]) . based on Ek 3> Aqcd, can only be viewed as a rough model 
calculation. The conclusion that Aag^H5{K) remains negligibly small should however still 
hold. The weak annihilation correction is Aag^^QaiK) = 0.003 at = 15 GeV^ for r = u 
and therefore entirely negligible, a consequence also of the small Wilson coefficients. 



4.2 B ^ K* 

In the case of the decay into a vector meson the relevant form factors are defined as 

(my = rriK*) 



{K*{k,e)\s^^h\B{p)) = -2i ^^^'^ 



TUB + my 



(33) 



[K*{k,e)\s-i^-i^h\B{p)) = 2mvAoiq^f—^q^ + (m^ + mv)A,{q 



2^ g • g 

rriB + my 



q 



(34) 
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It is convenient to treat the decay into longitudinally and transversely polarized 
vector mesons separately. Omitting terms proportional to 



rriB 



l^)h\B{p)) = 



2mv y/1 - {mv/EY mvimB + mv) 



A, 



(35) 



{Kl{k,e)\-sY{l - l,)h\B{p)) = '^-^e^^'-'el.ppk^ - {ms + my)A^e*l (36) 

In the large energy limit E ^ my, which we may use in the normalization of the 
power corrections, fl35|) and fl36l) simplify to [281129] 

(ir||(A;,£)|s7'^(l-75)6|5(p)) = -2A;Mo 

- 2V 

(iri(fc,£)|s7^(l-75)&|5(p)) = [^e^^P'^e\,v,k,^k■ve*l) (37) 

vrtB 

The case of a longitudinally polarized K* is very similar to the case of a pseudoscalar 
discussed in section 14.11 and we find 



Aa9,H5(i^ll) = ^A^^^^-/(") A,Ao(a^)a^ ^''^ 



where /y is the decay constant of K^^ 

For a K* with transverse polarization we obtain 

Aag^H5{K^) = — Ab 1/(g2) q'^ ^'^^'^"^^^^ + ^C^Qb) (39) 

where /_l is the decay constant of K]_. 

Numerically we have Aag^Hh{Kl) = 0.021 - 0.012i and Aag^n^iKl) = 0.008 - 0.006z 
at g2 = 15 GeV^ for our standard set of parameters. These corrections are of similar size 
as for the pseudoscalar kaon and they are likewise negligible. 

Finally, we quote the corrections from weak annihilation 

AW*I)^-(c..|)?g^ (40) 



A..uAi)^-(c,.|)?^(Q,-i) m 

where r = u, d refers to the spectator quark in the B meson. The corrections are 
tiny, at = 15 GeV^ we have AaQ^H6a{K*\) = Aag^H6a{Kl) = 0.003 for r = u, and 
Aag^HGaiK*^) = —0.001, Aa9^H6a(-^l) = —0.006 for r = d. This is again a consequence 
of the 1/g^ suppression at large values of g^, see also |30j. 
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5 Large vs. small recoil energy of the kaon 

The OPE for the correlator ( ITOl) . apphed to B K'^*h'^l^ , is vahd as long as the energy 
of the kaon in the 5-meson rest frame 

= < + (42) 

is small compared to This condition is certainly fulfilled in the vicinity of the 

endpoint, but even for as low as 15GeV^, just above the narrow-resonance region, 
Ek = 1.24 GeV is still fairly small in comparison to the hard scale. On the other 
hand, such a value of Ek is already larger than the QCD scale Aqcd and one could 
consider using the factorization methods applicable to the case of energetic kaons. For 

at 15 GeV^, or somewhat above, we have a transition region for the applicability of 
factorization methods for large Ek and the OPE for small Ek- It is this transition we 
wish to explore in the present section. 

QCD factorization (QCDF) at large kaon recoil requires Ek ^ Aqcd for arbitrary 
g^. The OPE method requires ^/(f ^ Ek, Aqcd- Both scenarios are consistent with 
the case 

V^? Ek:$> Aqcd (43) 

which can be realized, at least approximately, as we have seen above. In both, the QCDF 
and the OPE scenario, the leading contributions to the B — )■ K^*H^l~ amplitudes are 
given by short-distance quantities multiplying the standard B — t- fC*^*) form factors. 
In QCDF, the first corrections to these form-factor terms come from hard-spectator 
interactions, which have been computed in [TB]. Normalized to ag (called Cg in [15j ) 
these corrections may be written as AogY^ + Aag^l, where the indices refer to the 
notation of [15j. For definiteness we will consider the case of a pseudoscalar kaon in 
the following. Neglecting the small penguin coefficients C3, . . . , Cg and adapting the 
expressions in [15] to our notation, the second contribution reads {u = 1 — u) 

rrasCpmbfBfKQq f du (t>^{u) 



2N niBf+iq'^) J niBOJ — — ie 



^ 8Cs ems 



du (pKiu) 



U + M-V "^6 



H Cih{z,s)\2_ 



,,2 ^ ^ )\q'^^um%,+uq'' 

_2 rri A. ^ B ' 1 



(44) 



In the OPE limit, defined by treating g^ ~ and taking (H3ll . we have Aog^^ ~ l/m%- 
This term is therefore subleading with respect to (13T]) . consistent with the absence of a 
Cs term in Aag^HbiK)- The remaining term is given by 

a4Y) = -^C,Q,^¥^ f -Mu^) f'duMu)h (45) 



9,11+ 2N '^'mB/+(g2) J to 
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where {E = Ek) 



-#^1 + ^2E2 (Boiuml + uq') - B,{q')) (46) 



B,{s) = \ll-^^[ In ' + ZTT I (47) 



1/1 2 \ 1/2 1/1 2\ 1/2 

l./l mi \ . . 1 . / 1 mi^ 



z \4 um^ + uq^ / 2 \4 / 

Li{w±) = Lia ( "^"^ . 1 - TT + ln(l - w±) (in - — ^ ± zvr ) , w± = a;±, ?/± (50) 

The formulas from [IS] have been slightly rewritten here to make them more convenient 
for application in the high-g^ region where 4m^/g^ < 1. To make contact with the OPE 
result we eliminate q^ from (H6!) using g^ = — 2m bE and expand in E/mB- The 
terms proportional to l/E"^ and 1/E cancel and we find 

t,=/W + o(A) (51) 

with the function f{x) from ( l25i) . In the limit i^^i^- <^ ms the QCDF formula ( l45l) then 
reduces to the OPE result (I3T1) . We note that the hard-spectator correction (145|) . which 
is a leading-power contribution in QCDF, becomes a power-correction ~ 1/g^ ~ 
in the OPE regime. This behaviour is related to the form factor /_|_(g^), which scales as 
at small g^ and as m]/^ at large g^. 
We may use the preceding comparison to check the validity of the OPE at the rel- 
atively low values of g^ around 15 GeV^. Assuming that the large-energy limit for the 
kaon is a reasonable approximation, the QCDF expression (H5|) and the OPE result (13T1) 
differ only by the replacement 

/•I ^2 

/ rfM0^(M)tj, = -5.05 + 2.70i ^ /(x)^ = -5.87 + 3.55i (52) 
Jo 1 

where the numerical values are obtained at g^ = 15 GeV^ for our standard set of param- 
eters, neglecting the kaon mass for consistency. The real part of the OPE approximation 
is 16% larger in magnitude than the QCDF result, for the imaginary part the discrep- 
ancy is about 30%. Such differences are to be expected since the expansion leading from 
QCDF to the OPE is governed by Ek/ \f(p' ~ 0.3. Within this accuracy, the OPE for- 
mula still gives a very good estimate of the more complete QCDF result at g^ = 15 GeV^. 
For larger values of g^ the numerical difference between the two sides in ( |52ll becomes 
smaller and towards the endpoint the OPE is the more appropriate description. 
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The above exercise suggests that even at as low as 15 GeV^ the OPE is a vahd 
method to obtain the second-order correction (13T]) . The difference with the QCDF 
estimate is immaterial in view of the very small overall size of the effect. 

The transition from large to small recoil energy has already been considered in [22] for 
the B ^ K form factor ratio fx/ f+- Also in this and in analogous cases, hard-spectator 
corrections, which are leading-power effects at large recoil, become power suppressed for 
soft kaons. Similar observations hold for weak annihilation 1221. 



6 Toy-model analysis of duality violation 

The OPE defines a systematic framework to compute the correlator (fTOj) at high in 
QCD. In the Minkowski region > 0, there will be uncertainties in the OPE-based 
predictions that go beyond those due to neglected orders in as or in A/nib- Such effects 
are refered to as violations of quark-hadron duality. We investigate their importance 
first within a toy model for the charm loops in B — )■ Kl^l^ . In section 7 the model is 
generalized to be closer to the realistic case. 



6.1 Description of the model 

It is illuminating to consider the systematics of duality for the charm-loop contribution 
to i? — 7- Kl^l~ in the simplified context of the toy model introduced in [12]. This model 
assumes the existence of two leptons, li with a large mass mi and I2 with mass m2 = 0, 
and the effective weak Hamiltonian 

G - 

'Heff = [{yi)v-A {cc)v-A - (,yi)v^A {tt)v-A\ (53) 

All particles have standard strong and electromagnetic interactions. Then "Heff gives 
rise to a loop-induced process h — )■ I2 e^e~ via charm- and top-quark penguin diagrams 
with a GIM-like cancellation between them. The hadronic physics is fully contained in 
current correlators of the form 

nf = ^ I d'xe"^-^ (0|rj,^(x)j:(0)|0) = (g'^g^ - gV^) Il,{q^) (54) 

with = c'y^c. The decay amplitude can then be written as 

A{h ^ h e+e-) = --^ e,e' U{q') k^il - l^)h ei,e (55) 

where 11 = — II^ is the difference of the charm and top contributions. We take 
rrit > mi, and thus Imll comes only from the charm sector. The correlator 11 (g^) fulfills 
the dispersion relation 
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where the subtraction constant n(0) is fixed within the model and can be computed in 
perturbation theory. To leading order it reads 



n(o)^n.(o)-n,(o) = ^in!| (57) 

The form of fl55|) for the amplitude holds to lowest order in G and e^, but to all orders in 
the strong coupling. The decay li — )■ l2e^e~ in this model shares important similarities 
with B — )■ Ke^e~ , but the hadronic dynamics is simplified to the physics of quark-current 
correlators. The role of resonances and quark-hadron duality can thus be illustrated in 
a transparent way. 

From ( 155|) we obtain the differential decay rate (with s = q'^/mf) 

;?£(!i^ii£:£l= W(i_,).,i^2.)|C + A(,^)r (58) 

as IOStt^ ' ' 



where we defined 
and 

To lowest order we have 



C = 2n^n{0) (59) 
A(g2) = 271^ (n(g2) _ n(0)) (60) 



C = ln^ (61) 



and the partonic expression for A(g ) in one-loop approximation is 

r 1 I 2 arctan J—j ; a; > 1 

A,iq') = U^---A2 + x)y^\l^\\ (62) 
^ ^ 4 In -m- x<l 

where x = 4m^/g^. For typical values of the parameters (e.g. = 1.4 GeV, mi = 
mft = 4.8 GeV, rrit = 167 GeV) we have C = 4.78, whereas ReAg(g^) first rises from 
at = to 4/3 at = 4m^ and then drops again monotonically to the small negative 
value —0.07 at = m^. The coefficient C represents the short- distance contribution 
of the amplitude. It is real and larger (parametrically as well as numerically) than the 
quark- level charm contribution \Aq\. 

The decay rate (l58l) is proportional to 

|C + A|2 = c2 + 2CReA + |A|2 (63) 

As long as |A| is small compared to C, there is a clear hierarchy among the three terms 
on the r.h.s. of ( 163|) : The first, short-distance term dominates and the next term 
gives the correction to first-order in A, whereas the final term enters at second order. 
These features are qualitatively similar in the case of -B — )■ Kl^l~ . 
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6.2 Shifman model for charm correlator 



In order to investigate the systematics of duality violation in li — )■ l2e~^e~, we find it 
convenient to consider first a simple model for the quark-current correlator, which has 
been proposed in |311l32ll33] . In this model the correlator is represented as an infinite 
sum over resonances, which include finite width effects. In its original form it applies to 
massless quarks and we will correspondingly neglect the charm-quark mass in the present 
section. A detailed discussion of the model and its use in illustrating duality violation 
in the R ratio and similar quantities has been given in [3T]l32f33] . The model has also 
been used to study duality violation in r decays in [3^ . 
In Shifman's model the correlator A in fl60|l reads 

A(g^) = -f ^3^[V^(^ + l) + 7] (64) 
where ip{z) = dlnr{z)/dz is the digamma function and 

z = {-r~ ief-^l^ with r = f- (65) 

= 3 is the number of colours, A is a scale corresponding to the string tension in QCD 
and h = B/N = Tn/Mn is a (small) parameter related to the width-to-mass ratio of the 
resonances. 

The model expression for A in has the correct analytic behaviour (a cut for posi- 
tive g2 but no other singularities on the physical sheet) and it reproduces the asymptotic 
result of QCD in the limit of large g^, 

A(g2)^-^ln^^ (66) 

Using the identity 

ip{z + 1) + 7 = [ip^—z) + 7 — ivr] -j^ + [—71 cot nz + in]^ (67) 



the function in (iMll can be decomposed into two parts, A = Ai + A2, corresponding to 
the two brackets in fl67|) . A2 is an oscillating function of g^, exponentially suppressed 
for large g^. It represents the duality violating component of A. The function Ai is a 
monotonous function, which gives the OPE approximation to A. The real part of these 
two functions is shown in Fig. HI A plot of the imaginary part of A can be found in 
[32f33] . 

The duality violating part of ReA can be approximated as 

2\ ^ 1 A^vr 

ReA2(g ) = — — ^'^^ ^-^ + ^^J ~ — ^ exp{—2nbr) sin(27rr) (68) 

6 1 — o/vr 3 

if 27r6r ^ 1 and folnr/vr ^ 1. 
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Figure 4: Shifman model for charm loop: Re A(g^) as a function of /\^ for h = B/N 
1/6. The true function (oscillating curve) is compared with the OPE approximation. 



In the decay rate integrated over the high-g^ part of the spectrum, the duality vio- 
lating contribution enters proportional to 



ds {1 + 2s) {1 - sf Re A2 



Nn 



so 



ds (1 + 2s) (1 — sY exp(—2nbus) sin(27rus) 



N 



1 



— (1 + 2so)(l - So) - exp(-27rfooM) cos(27rsoM) + C -, — 
o u \u 



b 1 



where 



2 / 2 
q /m^ 



us 



(69) 



(70) 



In dnn]) we have used the approximation from fIBSl) . Typical values of the parameters are 

b = - u = 10 (71) 



1 

6 



The value of m = 10 corresponds for instance to = 2.3 GeV^ and ml = 23 GeV^. The 
quantity in f l69|) is shown in Fig. [5] as a function of sq. 

We comment on several important aspects of these results. 

• The duality violating component of ReA in fl68l) exhibits the characteristic oscillat- 
ing behaviour in r = with an exponential suppression governed by br. The 
analogous expression for the duality violating term in ImA, which has a cosine 
instead of the sine, has been given in 



Eq. (169 p displays the duality violating contribution from ReA to the partially in- 
tegrated decay rate. The integration over s extends from a suitably chosen lower 
limit So up to the end of the spectrum. The parameter sq should be large enough 
such that the OPE still remains reasonable at the corresponding value of q^. Using 
the approximation in fl68|) and expanding in the small quantities b and 1/m, we find 
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Figure 5: Shifman model for charm loop: Duality violating contribution J^|^ (l+2s)(l — 
s)^Re A2 as a function of sq with parameters h = 1/6 and u = 10. The full result (solid 
curve) is compared with the approximation given as the last term in eq. (l69l) (dashed 
curve). 



the explicit result written as the last expression in ( l69l) . This expression is an oscil- 
lating function of sq, with frequency 27iu, multiplied by an exponential suppression 
factor. The latter is active when the exponent is large, at least 27rbsou ~ 1. In 
addition, the entire term is further suppressed by one power of 1/ti = \^ /m\. This 
effect remains even if the exponential suppression is not fully developped. The 
power suppression arises because the oscillating contributions average out in the 
integral, except for a remainder ~ 1/m near the lower end of integration sq. 

• The duality violating term (IMj) is plotted in Fig. [S] for a semi-realistic choice of 
parameters, where mi = 4.8 GeV and = 2.3 GeV^. The smallest values of Sq 
shown, So ~ 0.05, correspond to ~ 1.2 GeV^ and 2TcbsoU ~ 0.5, which is already 
on the low side of the allowed range. From Fig. [5] we observe, first, that in the 
scenario considered here the simple approximation to ( !69l) agrees very well with 
the full result. Second, the numerical size of the duality violating term is about 
±0.02 for the relatively low values of sq around 0.1. The effect diminishes quickly 
for larger Sq due to the exponential suppression. The variation ±0.02 from (169|) 
amounts to ±6.5% when compared with the corresponding OPE expression 

[ rfs(l±2s)(l-s)2ReAi = -0.31 (72) 
Jo.i 

We next turn to a discussion of the |Ap term in fl63|) . where the situation is system- 
atically different from the case of Re A. The |Ap term receives a contribution from the 
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duality violating component given by 



ds{l + 2s){l- sf\A2\^ ^ J ds{l + 2s){l- sf y—J exp{-4nbus) 
■NttV {l + 2so){l-Sof 



4.5. exp(-4.6.o.) + O^^J (73) 

where in the second step an approximation similar to the one in ( 168|) has been used. 
There is still an exponential suppression, which makes the entire term negligible for suf- 
ficiently large BsqU. On the other hand, the power suppression with 1/u observed in fl69|) 
is softened into a behaviour as l/{bu). For small b (comparable to 1/u or smaller) the 
violation of duality may become large. This is in qualitative agreement with the discus- 
sion of duality violation for the squared correlator |n(g^)p given in [12]. The enhanced 
impact of duality violation is related to the absence of oscillations with alternating sign 
in the integrand of ( !73|) . We conclude that the | A p term in the integrated rate is partic- 
ularly susceptible to violations of quark-hadron duality, which may lead to substantial 
(positive) deviations from the OPE result, unless is large enough for a strong expo- 
nential suppression of the effect. However, the uncertainties in the |Ap term may be 
immaterial if this contribution is only a small correction to the dominant part in fl63|) . 

Let us finally illustrate the relative importance of the various contributions to the 
rate of li — )■ l2e'^e~ using the Shifman model for the charm loop with the numerical 
input defined above. In the limit — ?■ considered here the function C from fl^T]) is 
modified to 

C.U.^,l,l.,S2 (74) 
The relative contributions to the decay rate ( 1581) . (!63|) then read 



dsg{s)C^ = 13.60 (75) 



0.1 



[ ds g{s)2C Re Ai = -3.59 [±0.23] (76) 
Jo.i 

[ ds g{s) \Ai\^ = 1.25 [±0.10] (77) 

with g{s) = (1 ± 2s) (1 - sf. The central values in ([7SD and ([77]) are based on the 
OPE result for A, the square brackets indicate the impact of duality violation. The 
uncertainty in fl76]) gives the variation due to ReA2. It has the relative size of ±6.5%, 
which is reduced to ±2% in the sum of all contributions. The shift from duality violation 
in fl77]) is positive and essentially negligible in the present example. 

6.3 Model for charm correlator based on BES data 

In order to obtain a more realistic picture of the cc- spectrum, we are going to fit the 
BES data [5^361137] for the i?-ratio in the cc-region to a modified version of Shifman's 
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Figure 6: Simple fit to BES data for R vs. g^/GeV^. Riglit: Detailed View. 



model. The spectra of cc mesons can be accounted for by linear relations for the squared 



masses, = + M"^ 



' 



n 



1, 2, 3, . . ., similarly to the case of the light mesons 



The trajectory of the rt'Si charmonia, the J^^ = 1 states iIj{3097), tp{3686), '?/'(4040), 
'?/'(4415), . . ., for instance, follows this pattern. Starting form the third resonance (n = 3), 
these states can decay into open charm and have widths of the order of Aqqd- The 
first two are extremely narrow and may be described separately, but their properties 
are unimportant for duality violation, which is related to the infinite tower of high-n 
resonances. We therefore choose an ansatz where the sum over resonances begins at 
n = 3 rather than n = 1. A finite width is included in analogy to (!65|) and the variable 
is shifted by a constant into — 4m^. This leads to the following expression for the 
imaginary part of the correlator or, equivalently, the R ratio 



R — -Riip-ht — 



1 



3 (1 -6/7r)vr 



Im 'ip{3 + z) 



— Ami 



A2 



l-fc/vr 



(78) 



The individual resonances are located at = riA^ + 4m^ 



n 



3,4,5, 



in the limit 



6 — > 0. We observe that a rough description of the BES data [35] can already be obtained 
with this formula, where we find -Rught = 2.31 from the measured R-ratio below charm 
threshold, = 1.33 GeV, h ~ 0.082 and ~ 3.08 GeV^ This yields the result 
shown in Fig. El corresponding to a x^/d.o.f. ~ 2.5. We remark that our values for the 
fit parameters vtif. and A^ are in agreement with the results of [3H]- There is a second 
trajectory of 1 charmonia, the v?D\ states. Of these the first two resonances ?/'(3770) 
and '?/'(4160) are known. The first one is barely above threshold and still rather narrow. It 
may be considered separately, similar to ■?/'(3097) and ?/'(3686). Note also that ?/'(3770) is 
still below our default choice for the lower limit of the high-g^ region, > 15GeV^ The 
remaining resonances n^Di are rather close to the resonances {n + l)^Si for n > 2. For 
an approximate treatment it appears justified to subsume such a pair of close resonances 
under a single peak and keep the ansatz given in (!78|) . The accuracy of this description 
can be gauged by inspecting Fig. [61 In any case, the normalization of the second term 
of R in (1781) is fixed in the large-g^ limit by the free-quark result. 
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The fit could be refined by including the low-lying, narrow resonances in the fit ansatz. 
Any finite number of resonances does not change the asymptotic behaviour responsible 
for duality violation. Also the fit parameters are not much affected by such modifications. 
For example, including the ^'(3770) reduces the x^/d.o.f. — ?■ 1.7, while the parameters of 
the continuum ansatz remain almost unchanged, -Rught = 2.26, rric = 1.33 GeV, b ~ 0.078 
and ~ 3.08 GeV^. We will therefore be content with the simple representation given 
in f lTHj) above. 



7 Quark-hadron duality in B ^ Kl^l 
7.1 General considerations 

The analytic structure of the matrix element of the operator product in ( ITO|) can be 
infered from 

^'r(n) = I f d'xe''^-^ {K{k)\T f{x)H\Q)\B{p)) 



X 
X 



Pxo — Qq — 1^0 — 
Pxo + go - Po - «e 



{K{k)\H%0)\X){X\fmm) (79) 



In order to discuss the properties of the matrix element in (1791) we make the fol- 
lowing simplifications. First, we neglect the small penguin contributions in the weak 
Hamiltonian, that is we take 

= CiQI + (80) 

Second, in the electromagnetic current we retain only the charm-quark component, = 
QcC'j^c, and consider only contributions where the charm fields in H'^ are contracted with 
those in j^^. This neglects terms where the charm- anticharm pair from H"^ annihilates 
into gluons before connecting to the electromagnetic current. Such contributions are of 
higher order in and not essential for the problem we want to address. 

We consider the matrix element in fl79|) as a function of Qq, keeping the kaon energy 
fixed at a value of order 1 GeV. This can be achieved by injecting a spurion 4-momentum 
[39] r = (ro, 0, 0, 0) into the H'^ vertex. Then p + r = q + k and for po and ko fixed at their 
physical values the variable ro = go + ^0 — Po grows with go- The physical kinematics is 
recovered for r = 0. Under the simplifying assumptions specified above, the intermediate 
state X in the first sum in (1791) always contains a cc pair and a strange quark, in general 
together with other hadronic states, and thus pxo — ko ^ 2mc. The state X in the second 
sum contains a cc pair and a b quark, and pxo — Po ~ 2mc- The matrix element in (17^ 
is then seen to be an analytic function of go in the entire go-plane, except for two branch 
cuts at 2mc ^ Qq < 00 and at —00 < go ^ —2mc- If we would relax the simplifications 
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above and allowed for intermediate states X without cc pairs, the cuts would extend 
down to lower values of |go| on the real axis. 

The OPE of the matrix element in ( ff9l) can be justified for go on the imaginary axis, 
sufficiently far from the origin, that is at go = iQoE, for qoE ^ ^qcd- The OPE defined 
in this way in the Euclidean can then be analytically continued, term by term, from 
imaginary go onto the positive real axis, corresponding to the Minkowskian domain. 
Terms that are exponentially suppressed in Aqcd/Qoe for large positive go^; become 
oscillating functions of go in the Minkowskian case, that is for large positive go |39j . 
These oscillating terms are invisible at any finite order in the OPE and represent the 
duality violating contribution. 

7.2 Quantitative estimate of duality violation 

For a quantitative estimate of duality violation we have to resort to a model of the 
hadronic correlator in ( |79ll . To this end we write the Hamiltonian in (IHOl) as 

= a2 isb)v^Aicc)v-A (81) 

and assume a factorization of the currents, that is we neglect interactions between cc and 
the B ^ K system. The coefficient 02 is then treated as a phenomenological parameter. 
With these simplifications the correlator in (1791) reduces to 



1 Rtt^ 

«) = ^a2((^-&)y-A)^n,(g2) (82) 

where we omitted the longitudinal component ~ g'^. lie is the current correlator defined 
in dSlD- The charm loop in (1521) contributes to the coefficient Og in the amplitude of 
B Kl^l~ a term 

Aag = a2d, d = -— {U,{q^) - n,(0)) (83) 



In the model of section 16.31 we have 

4 



d=-\ ^ \, [^(^ + 3) - ^(^0 + 3)] (84) 



g^ — Ami 



where 

z = {-r - itY-''\ r='^^f^ = u{s~s,), u = m],/X' (85) 
and Zo = z{q^ = 0). For the parameters we use the following values 

= 3.08 GeV^ = 1.33 GeV, b = 0.082 (86) 



We will not employ the model (1831) to describe the entire charm-loop contribution, but 
only to estimate its duality violating component. The remainder is more reliably obtained 
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by the OPE itself. In order to extract the term that represents duahty violation within 
the model, we decompose 



tpiz + 3) - ^{zq + 3) = [ip{~z - 2) - 4j{zo + 3) - ivr]^ + [-vr cot nz + in]^ (87) 

This decomposition is useful for > 15GeV^, when the second term starts being ex- 
ponentially suppressed and gives the duality violating contribution. In correspondence 
with ( 1871) we have d = di + d2 and the duality violating term is 

1 r . , 871 

——— \ —7r cot TTZ + zrrl^ ~ 

3 1 - Vtt ^ 3 



d2 = —-- TTT '^ot + ^^]2 ~ — ~ exp(— 27r6r) (sin 27rr — i cos2nr) (88) 



When integrating the jagp part of the B — )■ Kl^l^ rate over the high-g^ region, from 
a lower limit = Sq^^ to the end of the spectrum, the relative size of the duality 
violating effect is then given by 

2a j::dsXf{s)flis)Red2 

J^DV,l — Z Tin (.oyj 

a, CdsXf{s)fUs) 

to first order in the charm-loop contribution. Here we have assumed 02 to be real. Using 
the approximation in flHHl) and proceeding as in fl69|) we find 

ID |/8«2 a5/^(so)/|(so) 1 / o /, r ^^ rnn^ 
\Rdv,i \ ^ 77 0/0 ^ -exp(-27r6M(go - sj) (90) 



u 



This formula gives an excellent approximation of the full result based on ( 189|) . We note 
the (mild) exponential suppression and the power suppression by l/w = X^ /m\. For a2 
we take the value 02 = 0.3, which is large enough to reproduce the measured B — )■ Kifj 
branching fraction within the factorization ansatz. We recall that ag 4. Then, for 
0.5 < So < 0.6, that is for q^ in the vicinity of 15 GeV^, the relative correction |-RDy,i| is 
below 3%. In the rate for B — )■ Kl~^l~, |aioP is added to jagp, which roughly doubles 
the result. The net effect of the uncertainty from (1901) for the rate integrated over the 
high-g^ region is therefore only about 1.5%. 

As discussed in section [6l the second order effect in ^2 is qualitatively different. It 
has no cancellations due to oscillating terms, giving a positive correction, and its impact 
increases with decreasing b. If local duality is at least roughly fulfilled, as is the case for 
high enough > 15GeV^, the duality violation from |(i2p is still suppressed, being of 
second order in the small quantity a2/ag. The relative size of this component is 



^ _4 LTdsXf{s)fl{s)\d2\^ 
J^DV,2 — —9 — 

al f:;;dsxf{s)fi{s) 



3«W CdsXfis)l 



(91) 
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where the second term uses the approximation in ( 188|) . The approximate form and the 
full result in (19T|) agree reasonably well. With our set of parameters we find Rdv,2 = 0.015 
at gQ = 15 GeV^, about half the size of |-RDy,i|- 

The phenomenological value 02 = 0.3 used above is larger than the perturbative result 
for the coefficient 02. Thus it effectively absorbs the (partly unknown) effects from fac- 
torizable and non-factorizable corrections. We should emphasize here that the analytic 
structure of the latter is, in general, more complicated than it is implied by the approxi- 
mation to (/C^) in f l82|) . Globally using the larger value of 02 in our numerical estimate of 
duality violation thus corresponds to the pessimistic scenario where the oscillating terms 
from the non-factorizable corrections to the charm-loop are added coherently. In reality, 
we expect that at least some destructive interference between the various contributions 
appears, and the amount of duality violation should even be smaller than our estimate. 

Our model for the charm-loop in ( 182|) is similar to the ansatz originally proposed in 
|4U] and used since then in many phenomenological studies. However, our motivation 
for considering this model is essentially different. In contrast to ^D] it is not our goal 
to model the hadronic effects on the spectrum point by point in the distribution. 
More relevant than the detailed shape of the spectrum is the rate integrated over the 
entire high-g^ region, which is best described in a model-independent way by the OPE, 
as mentioned above. We rather employ the model to get an indication of the duality 
violating effects, which are not captured in an OPE calculation. For this purpose the 
Shifman model for the current correlator is adapted to the charm-quark case He, with a 
choice of parameters consistent with the basic features of the most recent experimental 
data from BES. Note that the Kriiger-Sehgal ansatz for 11^ |1Q], consisting of a spectral 
function with a finite number of resonances and a flat continuum for large g^, contains no 
information on duality violation. The model we use to estimate duality violation is very 
simple and involves many assumptions. Still we expect it to indicate the systematics and 
the typical size of the effect, which presumably is closely connected with the cc resonance 
structure in the charm-loop contribution. 

8 Comments on the literature 

The high-g^ region of 5 — )■ K*l^l~ has also been analyzed in the framework of an OPE 
in [18] and a recent application was presented in [H]. Our approach differs from the 
analysis of [18] in several respects. We go beyond the work of [18] by addressing in detail 
the issue of duality violation and the basis of the OPE formalism. In addition, we extend 
the OPE to include second-order power corrections, which we estimate quantitatively. 
A list of new results is given in the Conclusions. Here we would like to comment further 
on two conceptual differences between [18] and our formulation. This concerns, first, the 
construction of operators in the OPE and, second, the treatment of charm quarks. 

The operators in our approach are built from 6-quark fields in full QCD rather than 
using HQET. This is convenient because the operator basis is simpler and the OPE 
becomes particularly transparent. Another advantage is that the matrix elements of the 
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leading operators are given by the usual form factors in full QCD. Unlike in a HQET 
framework, not all dependence on nif, is made explicit, but this is not essential from a 
practical point of view and still allows the consistent inclusion of power corrections to 
a given order in the expansion. This approach has been used for instance in the OPE 
for inclusive B decays applied to the computation of the lifetime difference of Bg mesons 
[^2|l3] . These investigations included power corrections ^2] as well as corrections of 
0{as) [l3]. By contrast, in |18] a matching onto HQET operators is performed from 
the start. This leads to a proliferation of operators, whose matrix elements are not 
given by the usual form factors. In fact, to reach a simplification of the resulting HQET 
expressions, the authors of [I8] partly undo the matching , from HQET back to full QCD 
expressions. We prefer to employ the fuU-QCD formulation throughout, in view of the 
advantages mentioned above. 

We next turn to the second point, the treatment of charm. The authors of [18] 
perform an explicit expansion in rric/mf,, which corresponds to assuming the hierarchy 
rric <^ rri!,. This implies that the charm quark is not integrated out at the scale mt, ~ ^/(f 
and continues to be an active field below this scale. Operators with charm-quark fields 
are then present in the OPE. 

We will argue that it is conceptually simpler to integrate out charm immediately at 
the scale rrib, assuming a hierarchy nib rric ^qcd, and that this can be done without 
any loss in accuracy. Charm-quark effects are then entirely contained in the coefficients 
of the operators, as is apparent from the formulation of the OPE given in section 3. 

To illustrate the point we consider the following example. Let us take the scenario 
where ^ nib. In this case the operator coefficients should be evaluated with ni^ = 0. 
On the other hand, additional operators involving charm would have to be included in the 
OPE. With explicit cc fields, such operators arise at dimension 6 or higher. For instance, 
radiating the virtual photon from the two charm lines in the effective Hamiltonian (IHOll 
and leaving all four quark lines open, gives the dimension-6 operator 

i5ij5kiCi + 5ii5kjC2) (92) 

in analogy to fl26|) . This operator has B — K^*^ matrix elements contributing to the de- 
cay amplitude. Assuming next that ^ ^qcd, the contribution fl92|) can be simplified 
by integrating out the charm fields in a further step. Contracting the charm lines into a 
loop, to which a gluon field is attached, induces below nic the dimension-6 operator 

^^6d = tCiQc In 4 ^'^''^^ - 76)9D^G,f,b (93) 
o q 

The coefficient of this operator reproduces the logarithmic mc-dependence at this order in 
the OPE. A similar discussion has been presented in [TH], illustrating how the logarithmic 
term In in the coefficient of the leading dimension-3 operator (122|) is recovered in an 
effective theory with four-quark operators of the type (s6)(cc). Whereas the latter effect 
vanishes for mc — )■ 0, we note that the coefficient in (p3l) has a logarithmic divergence 



Q7''(l - 75)6j Sfc7^(l - 75)q - {/i A} 
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in this limit. This is of no consequence if rric ^ ^qcd, where rric still represents a hard 
scale. In fact, for higher- dimensional operators generated by ( 1920 . that is operators with 
more factors of the gluon field and its derivatives, the coefficients will scale as inverse 
powers of rric. A similar situation exists for inclusive semi-leptonic b c decays, where 
four-quark operators with charm, analogous to (1921) . also appear at third order in the 
OPE for rric ^ ^6- The effect of the corresponding charm loops has been refered to as 
'intrinsic charm' in pi|45|46] . where the issue was discussed in great detail. The above 
consideration indicates how the nonperturbative B — K* matrix element of operators 
such as ( 192|) can be treated in an expansion in Aqcd/^c- However, as demonstrated 
in (iHl45ll46j for semi-leptonic b ^ c decays, the same effects are also described in a 
framework where charm is integrated out at the scale. No active charm fields need 
to be considered in this case and (!92|) is absent from the OPE. A difference between the 
OPE with or without active charm is that the framework with charm fields and with a 
strong hierachy rric ^ assumed, would offer the possibility of efficiently resumming 
logarithmic terms lnmf,/mc. Since such logarithms are not very large and in view of 
the additional power suppression of such terms, such resummations appear not to be 
necessary in practice. We also stress that effects such as f l92|) . irrespective of their 
detailed treatment, are suppressed at least as l/mf. Small differences in the method of 
their calculation are therefore hardly relevant. We thus conclude that integrating out 
charm at the scale mi, is entirely justified. The OPE is then constructed in a single step 
and with a simpler operator basis. For these reasons the approach appears preferable 
and we have adopted it here. 

9 Conclusions 

The amplitude for B — > Ml^l~ , M = K, K*, . . ., has the general form given in ([8]). It 
contains the component 



= C, {M\sr{l - ln)b\B) + Cj -^qx {M\sa'^{l + j,)b\B) + (M|/C^(g)|5) (94) 



which receives a nonlocal, hadronic contribution (/C^) = (M|/C^|i?) from the matrix 
element of the nonleptonic weak Hamiltonian H'^ in addition to the semileptonic form 
factor terms (~ Cg, C7). Although the hadronic part (/C^) is relatively small numerically 
outside the narrow-resonance region (~ 10% of Ag), it needs to be reliably computed in 
order to achieve very accurate predictions. We have presented a detailed study of the 
hadronic contribution in the region of large dilepton invariant mass > 15GeV^ based 
on an operator product expansion in inverse powers of the hard scale ^/(f. Working 
with 6-quark fields in full QCD and factorizing the dependence on rrib, \f(f and rric into 
the coefficient functions, we obtain the following results: 




(95) 
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• To leading order in the OPE and to all orders in a^, (/C^) is expressed in terms 
of the standard form factors parametrizing the matrix elements (s7'^(l — 75)6) and 
(scr^^(l + 75)6) of dimension-3 operators, up to coefficients calculable in perturba- 
tion theory. To lowest order in as only (57^(1 — 75)&) is present. 

• In the chiral limit (m^ = 0), the first power corrections appear only at second order 
(~ l/o'^) and are governed by dimension-5 operators with a gluon field of the form 
sgGVh. The corrections are computed explicitly, using the hmit Ek ^ ^qcd for 
the hadronic matrix elements, and are shown to be smaller than 1% for Ag. 

• For nis 7^ the dimension-4 operators msS7^(l + 75)6 and msSa'^^il — 75)6, with 
right-handed strange quarks, can arise in the OPE. Because they are absent at order 
ftg, their contribution will be suppressed to the negligible level of a^ms/mfe ~ 0.5%. 
Besides, no new form factors will be introduced by these operators. 

• Within the OPE framework the effect of weak annihilation is a natural ingredient, 
which we have briefly discussed, mainly for illustration. Since in addition to being a 
third-order power correction it comes with small Wilson coefficients, its numerical 
impact of a few permille is entirely negligible. 

• We have clarified the relationship between the OPE at high and QCD fac- 
torization at low by showing that both descriptions yield consistent results at 
intermediate values of q^ ~ 15 GeV^. 

• A relevant topic at high q^ is the issue of quark-hadron duality, which is closely 
related to the existence of an OPE. We have defined the OPE with the help of a 
spurion momentum that allows for an independent scaling of q^ at fixed itlb and 
kaon energy. This allows one to clarify the analytic structure of the matrix ele- 
ments of the correlator (p5|) . We then employed a model based on an infinite series 
of charm resonances to estimate quantitatively the amount of duality violation, 
resulting in about ±2% for the rate integrated over the high-g^ region. An impor- 
tant aspect is the different sensitivity to duality violation of the contributions to 
the rate linear or quadratic in (/C^). The quadratic term is more vulnerable to 
duality violation, but numerically suppressed as a term of second order in the small 
charm-loop contribution. The systematics of duality have further been studied in 
a toy model for the rare decays, in which the factorization of the charm loop is 
exact. 

The main conclusion is that the high-g^ region of i? — K^*H~^l~ is theoretically 
under excellent control. Decay rates and distributions are perturbatively calculable up 
to the nonperturbative effects accounted for by the standard form factors in full QCD. 
Further nonperturbative corrections are strongly suppressed and negligible within the 
accuracy of a few percent. The existence of the OPE implies that at high the theory 
ofB^ K^*H+l- has an even more solid basis than at low q^. An example for applications 
is the combined analysis of S — )• Kl^l~ and B — )■ Kuu [22P7] . Here the form factor 
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dependence can be essentially eliminated, which leads to precision observables sensitive 
to new physics effects [22]. Beyond B — )■ K^*H^l~, the results of our analysis apply to 
the high-g^ region of exclusive rare decays with similar final states such as -B — ?■ Knl^l^ , 
Bg —7- (pl^l^ and, with appropriate modifications, to i? — > pl^l^ or i? — t- nl^l^ . 

A Complete basis of operators through dimension 4 

We show that the OPE of /C^ in ( ITOl) can be expressed in terms of 

^3,1 = {a"" - ^) ^7.(1 - 75)& , Ol, = ^g, sa'^^il + 75)6 (96) 

at the level of dimension-3 and dimension-4 operators in full QCD, in the chiral limit 
rris = 0. For rris 7^ 0, ^ rrib, right-handed strange quarks can also appear, together 
with a suppression factor of rris/mi,. This leads to the two additional operators O^i and 
C^4 2 f|T7|) and f|T8l) . which can be counted as operators of dimension 4. The matrix 
elements of O3 1, (^32; and O42 between B and K^*^ are all given by the standard 
B — )■ K^*^ form factors. 

Together this implies that at leading (dimension 3) and next-to-leading order (di- 
mension 4) in the OPE, and to all orders in a^, standard form factors are the only 
hadronic matrix elements required. Corrections to these terms arise only at second or- 
der (~ A^/q^) in the OPE through operators of dimension 5. In other words, there 
are no first-order power corrections (~ A/rrif,) in the OPE of A^^, except for the purely 
kinematical dependence on and ml in the coefficient functions. These functions con- 
tain first-order terms such as (m^ — q'^)/ml, which however are calculable and do not 
introduce unknown hadronic matrix elements. 

We demonstrate the completeness of the basis C3 1, £^3 2, C4 1 and C4 2 for opera- 
tors of dimension 3 and 4 by enumerating all possibilities consistent with the relevant 
symmetries, and using the equations of motion. 

We first assume rris = 0. Operators of dimension 3 built from sl and b have the form 

{sLTbr, r = 1,7^,0^"^ (97) 

where Lorentz indices can be contracted with the metric g, the e-tensor and factors of q 
to yield a 4-vector with index /x. Useful relations are 

q^'isLTb) =rdf'{sLTb) (98) 

which holds for B{p) — > K^*\k) matrix elements, and 

d^isLTb) = SLv'D^b + SLTD^b (99) 
or similar identities. Exactly three structures can be formed: 

SLl^b, quSLO^'b, q^Lb (100) 
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This can be seen as follows. For F = 1 in ( 197|) the only possibility is the third operator 
in ( 11001) ■ For F = 7° we obtain the first operator in ( llOOp and 

q^quSiYh = qf'isLi + Sii pb) = mbq^'sLb (101) 

which is equivalent to the third term in fllOOp . For F = cr"'^ we can write the second 
operator in (llOOp . In addition, we could form 

e>"'^''qo.SLcyxub ~ qaSicr'^^l^b = qaSL^^'^b (102) 

which leads back to the same structure. This exhausts the possibilities and proves the 
completeness of the basis in fllOOp at the dimension-3 level. 

We next consider operators of dimension 4. These have, in addition to the fields in 
( l97|) . one covariant derivative acting on the strange quark. Their general form is 

{sL^^Tbr (103) 

Because of 

SL^^VDxb = dxisL^Tb) - sl^^aF6 = -iqx{sL^Tb) + dim 5 (104) 

terms with extra covariant derivatives acting on b do not lead to independent dimension-4 
operators and may be disregarded. 

For F = 1 in f llOSp we can write the operator Using the identity in (fT6ll this 

operator can be expressed in terms of the basis in f llOOp : 

= '-m,SLrb - ^q^SLcr^'^b - '-q^Ts^b (105) 
The other possible structure can be reduced as 

q'^q.SL^'b = '-{ml - q^)q^SLb (106) 
For F = 7" the dimension-4 operator has the form 

SLl5''-f"b (107) 
We list the possible contractions of f ll07p into a 4-vector together with the reduction to 
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the basis in fllOOl) . neglecting operators of dimension 5. 
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qi'sL pb = (108) 
q^^q^q^SL^^'^'^h = '-mh{ml - q')q>'sLb (109) 

= ^d'-d^isL^b) - -SLi^^D'-D.b = -JL_J_SLl^b (110) 

qasj^f'-f'^b = slIdH pb = mbSLl5^'b (111) 

^''"'qaSL'5xlub = qaSL^Sxil'rr - 9^^l" ' Q""^!^ + g'^^l^h = 

= -q^SL^S^^b + q^SL^x^b (112) 

The r.h.s. of fillip is reduced to the basic operators via (11051) . The terms on the r.h.s. 
of (I112p lead back to the expressions fillip and flllOp . 
For r = cr"^ the dimension-4 operator becomes 

sl^V"'^6 (113) 
We list again the possible contractions and their reduction to the operators in fllOOp . 

sil^xcr^^b = isL^^xil^Y - g^^)b = -isL^% (114) 

qxqpSLl5^a^^b = im,qxSLl5^^''b - iqxq^SL^^^b (115) 

£'^^"%L^Acra/3& ~ SL^^xa^^b (116) 

The r.h.s. of ffTTip is given by flTOSj) . The r.h.s. of flTTSj) is equivalent to ffTTOj) and ffT06|) . 
and flTTej) reduces to ffTTip . 

This completes the proof that all operators of dimension 3 and 4, fp7|) and fll03p . 
can be expressed in terms of f llOOp for = 0. Taking current conservation into account 
leaves us with Cg^^ and C3 2 in (jHS])- Similar arguments hold if sl is replaced by sr. 
Since right-handed strange quarks come with a factor m^, the two further operators O'^^ 
and C4 2 in f ITTP and f|T8|) are obtained for 7^ 0. 

Similarly to the derivations above it can be shown that dimension-5 operators of the 
form slD i^Vb, contracted to a Lorentz vector with index /x, can always be reduced to 
linear combinations of f llOOp and operators containing a factor of the gluon field strength 
Gap- The former terms correspond to dimension-3 operators with a purely kinematic 
power suppression. Genuine operators of dimension 5 therefore have the form [siGTb)^, 
as quoted in f !T9|) . 
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rriK 


h 




f± [50] 


0.496 


0.16 


0.218 


0.185 


rriB 


TUB* 


/b 


Ab 


5.28 


5.41 


0.2 


0.350 


rhb 


rfic 




rriK* 


4.2 


1.3 


0.225 


0.894 



Table 1: Input parameters in GeV. f± = f±{l GeV). 



B Numerical input 

In this appendix we collect input we have used in numerical calculations. The numbers 
quoted are our central values. The form factors have an uncertainty of roughly ±15%. 
The B — )■ K* form factors are parametrized as 



1.364 0.990 
gV(5.28GeV)2 ~ 1 - g7(36.78 GeV^ 



Ml') = . „..;.oo^..n. - . .o^,..2^ (117) 



0.923 0.511 
1 - gV(5.32GeV)2 " l-g2/(49.40 GeV 

The B ^ K form factor is parametrized as [22 



with the default choice /+(0) = 0.304, ao=1.6 and bi = Bq. 
Further parameters are summarized in Table [H 
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